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1. Introduction 

The possibility of realizing sub shot noise regime by exploiting multimode spatial correlation 
at the quantum level |[Tl|2l, sometimes called spatial entanglement, has been experimentally 
demonstrated by using traveling wave parametric amplifier and CCD array detectors both at 
the single photon level fS) and at hight photon flux ||4l|5]|6ll2][8l, and by exploiting four wave 
mixing in rubidium vapors ifTOl . Very recently, following the proposal of |f9l our group showed 
that it can find a natural application to the quantum imaging of weak absorbing object beyond 
the shot-noise-level ||6|(standard quantum limit). All these findings indicate that multimode 
spatial correlations and their detection is a mature field that can lead to other interesting ap- 
plications. In particular, in this framework we developed the capability of an almost perfect 
spatial selection of modes that are pairwise correlated at the quantum level in parametric down 
conversion, and the optimization of the noise reduction below the shot-noise-level. 

On the other side, quantum correlations of twin beams generated by Parametric Down Con- 
version (PDC), is a well recognized tool for calibrating single photon detectors, competing 
with the traditional ones of the optical metrology nT1[ni[T3l[T4l[T5l[T6l[T7l[T8ll20ll2T1l221l23]|. 
Extending this technique to higher photon fluxes, for calibrating analog detectors, may have 
great importance in metrology ll24l l26l l27l l25l [191 l28l . So far, the main problem to achieve this 
goal has been the difficulty of an accurate spatial selection of the correlated modes among the 
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Fig. 1. Scheme for spatial correlations detection in PDC. The PDC emission from the non- 
linear crystal is detected in the far field, reached by a optical / — / configuration. The 
pump transverse size determines an uncertainty in the photon propagation direction. The 
speckled image showed, has been obtained experimentally in very high gain regime. For 
this reason, the spatial fluctuations are very strong, and the coherence area is roughly 
represented by the typical size of the speckles. 



spatially broadband emission of PDC. 

In the work presented here, the know-out in the detection of spatial correlation in PDC is 
fruitfully applied to the absolute calibration of CCD cameras. The method is based on the 
measurement of the degree of correlation between symmetric areas belonging to the twin 
beams, which in principle depends only by the transmission and detection efficiency. The 
achieved statistical uncertainty indicates the effectiveness of the methods for metrology ap- 
plications. 

2. Multimode spatial correlation in PDC 

The state produced by spontaneous PDC, in the approximation of plane wave pump field of 
frequency cOp and propagating in the z direction, presents perfect transverse momentum phase- 
matching. Thus, it can be expressed as a tensor product of bipartite states univocally identified 
by frequency C0p/2± co and transverse momentum ±q, i.e. \^') = 0q <„ \\ff{q,(o)). Since we 
are mainly interested to the frequencies near to the degeneracy (CO ^ 0), the state of the single 
bipartite transverse mode reduces to 



IWq)>=Ec"q(«)l«)'.ql«)^-.-q' 



(1) 



where the coefficients Cq{n) °^ \/("q)"/("q + 1)"^' are related to the mean number of photon 
in the mode q assumed to be the same for all the modes, i.e. («q) = jj., and the subscripts s and 
/ indicated signal and idler fields. 

The two-mode state in Eq. ([T]l is entangled in the number of photons for each pair of modes 
±q, whereas the statistics of the single mode is thermal with mean value («,.q) = (n.s.-q) = jti 
and variance (5^n,-.q) = {d^iis-q) = /x(/X + 1). Now, we focus on the far field region, obtained 
as the focal plane of a thin lens of focal length / in a / — / configuration (Fig. |2]l. Here, any 
transverse mode q is associated with a single position x in the detection (focal) plane according 



to the geometric transformation {2cf/cOp)q -^ x, with c the speed of hght. Therefore, a perfect 
correlation appears in the photon number «, x and «s._x registered by two detectors placed in 
two symmetrical positions x and — x, where the center of symmetry (CS) is basically the pump- 
detection plane interception (neglecting the walk off). 

In real experiments the pump laser is not a plane wave, rather it can be reasonably represented 
by a gaussian distribution with spatial waist Wp. This induces an uncertainty in the relative prop- 
agation directions of the twin photons of the order of the angular bandwidth of the pump. This 
uncertainty is the coherence area of the process, roughly corresponding to the transverse size 
of the mode in the far field (see Fig. |2]i. The number of photons collected in symmetrical por- 
tions of the far-field zone are perfectly quantum correlated only when the detection areas si/det 
are broader than a coherence area, whose size £/coh is of the order of [{2ncf)/{(0pWp)]'^ illlD 
(only for very large parametric gain it deviates from this behavior ||9l|4l|29l[30l). Therefore, 
let us consider to collect photons over two perfectly symmetrical and correlated areas £/det,s 
and £/jet,i with detection efficiency rjs and 77, belonging to the signal and idler beams respec- 
tively, and containing many transverse spatial modes .//i^spatiai = ^det^jl s^mh (j = ^, 0- We also 
consider a situation in which the detection time 3det is much larger than the coherence time 
^cohof the process, thus the number of temporal mode is large, ^t = 3'detl S'coh ^ 1- Since 
the modes in the single region are independent, the statistics is multithermal with mean value 
(A^;) = ^un'^jV^ (with Jitot = -^t-^spatiai) and variance 

{8^Nj) ^ (Nj) {l+^)^{Nj)ll + ^\^ J^„„i]j^ (1 + 77,-M) , (2) 

with i = s,i and S the excess noise, usually defined as the fluctuations that exceed the shot 
noise level (SNL). The SNL, or standard quantum limit represents the level of noise associated 
to the detection of the coherent light, i.e. {5^Nj)snl = {Nj). In this theoretical description, 
the excess noise is only related to the intrinsic thermal statistic of the single beam of PDC 
# = {Nj)/^tot- However, we will discuss in the following that experimental imperfections 
give the major contribution to the excess noise in our setup. 

The covariance between the signal and idler numbers of photon is 

{dNiSNs) = ^,„rTJ,T7;^(l + ^). (3) 

The amount of correlation between the signal and idler fields is usually expressed in terms of 
the noise reduction factor a defined as the fluctuation of the difference A^_ = Ns — Ni between 
the photons number normalized to the corresponding level of shot noise; 

(5^N\ ri^ fl \ n^ f {Ns+Ni)\ 



{Ni+N) ■'"' 277+ V2 "; '" 4rj 

where 77+ = (t], + T7,)/2 and ?]_ = r\s — ^i- It as been evaluated by introducing the Eq.s (O and 
© in the expression of the fluctuation {5-N-) = (S^N,) + {5^Ni) - 2(8Ni8N,) 

In the case of perfect balanced losses rjs = '?! = '? one get that a = 1 — 77 only depending on 
the quantum efficiency. Therefore, in an ideal case in which 7j ^^ 1, CJ approaches zero. On the 
other side, for classical states of light the degree of correlation is bounded by a > 1, where the 
lowest limit is reached for coherent beams, (7=1. 

According to Eq. (HI, the absolute estimation of the quantum efficiency by measuring the 
noise reduction factor, can be achieved just when the excess noise disappears, that is realized in 
the case of balanced losses. The balancing of the two channels can be performed physically by 
adding proper absorbing filters in the optical paths. Anyway, a more convenient approach is to 



compensate a posteriori, for instance multiplying the values of A^,- for a factor a = (Ns) / (Nj) = 
Tjs/Tjj. It corresponds to evaluate the redefined noise reduction factor Ga (instead of the one in 
Eq. dUl): 

This relation shows that the quantum efficiency rjs can be evaluated measuring Ga and the ratio 
a = (Ns) I (Ni), without the need of a physical balancing the two optical paths. 

As a final remark, we observe that the result in Eq. ^ relays on the assumption that each 
spatial modes collected by the region s^det.i finds its correlated in the region s^det.s^ and vicev- 
ersa. Otherwise, the presence of uncorrelated modes in the two region would not provide a 
complete cancelation of the excess noise by the subtraction, leading to underestimate the quan- 
tum efficiency ||3T1 . Therefore, experimental control on the detection of the spatial modes by 
means of precise positioning and sizing of the regions and accurate determination of the CS is 
fundamental for the accuracy of the estimation of the quantum efficiency. 

3. The experimental procedure 

In our setup, a type II BBO non-linear crystal {l — l mm) is pumped by the third harmonic (355 
nm) of a Q-switched Nd:Yag laser. The pulses have a duration of ^ = 5 ns with a repetition 
rate of 10 Hz and a maximum energy, at the selected wavelength, of about 200 mJ. The pump 
beam crosses a spatial filter (a lens with a focal length of 50 cm and a diamond pin-hole, 250 
/xm of diameter), in order to eliminate the non-gaussian components and then it is collimated 
by a system of lenses to a diameter of Wp = 1 .25 mm. After the crystal, the pump is stopped by 
a couple of UV mirrors, transparent to the visible (~ 98% transmission at 710 nm), and by a 
low frequency-pass filter (~ 95% transmission at 710 nm). The down converted beams (signal 
and idler) are separated in polarization by two polarizers (97% transmission) and finally the far 
field is imaged by a CCD camera. We used a 1340X400 CCD array, Princeton Pixis:400BR 
(pixel size of 20 /im), with high quantum efficiency (around 80%) and low noise in the read 
out process (4 electrons/pixel). The CCD exposure time is set by a mechanical shutter to 90 ms, 
thus each image acquired corresponds to the PDC emission generated by a single shot of the 
laser The far field is observed at the focal plane of the lens with 10 cm focus in a / — / optical 
configuration. 

For reasons related to the visibility of the correlation, and in order to reduce the contribution 
of the read noise of the CCD, it is convenient to perform an hardware binning of the physical 
pixels. It consists in grouping the physical pixels in squared blocks, each of them being pro- 
cessed by the CCD electronics as single "superpixel". Depending on the measurement, the size 
of the superpixel can be set accordingly. Typically we choose it of the same order, or larger than 

S^coh- 

The expected number of temporal modes ^t = ^pl -S^coh detected in one image is 5 • 10^, 
considering the coherence time 3^coh of PDC around one picosecond. The number of spatial 
modes ^spatial = -s^detj/s^coh ij = ■*, depends only on the size of the detection areas, since 
£/coh ~ 120 X 120(/im)^ is fixed by the pump transverse size. The level of excess noise due 
to the thermal statistics in the single beam is also fixed, since we keep fixed (aside unwanted 
fluctuation pulse-to pulse) the power of the laser. The total number of modes ^tot turn out to 
be compatible with the level of excess noise S" = (Ns) j ^tot '--^ 0.1 — 0.2. It can be measured 
by performing spatial statistics as described in Sec. (13. 11 1. 

For an accurate estimation of quantum efficiency the following steps should be performed: 

a) Determination of the center of symmetry (CS) 

positioning of the correlated areas and determination of the center of symmetry of the 



spatial correlations within sub-coherence-area uncertainty, according to the experimental 
procedure presented in Subsection ( 13. Il l 

b) Determination of the minimum size of s^j^f j 

The size of the detection areas must satisfy the condition ./^spatial = ^det.jl ^coh ^ 1> 
for the purposes of the unbiased estimation of O. This can be achieved following the 
procedure sketched in Subsection ( I3.2l i. 

c) Analysis of experimental contributions to the excess noise 

Evaluation of noise coming from experimental imperfections, such as instability of the 
laser pulse-to-pulse energy and the background due to straylight and electronic noise of 
the CCD. Eq.s ^ and (|7]l are modified in order to account for these noise contributions. 
Detailed discussion on this item can be found in Subsection ( I3.3l l. 

d) Estimation ofrjj and of its statistical uncertainty according to Eq. Q 

a and Oa are estimated experimentally over a set of ^ images according to formula 

a = ^M^ (6) 

and formula 

_ £[(M.(fc) - aNi{k)y-]-Emk) - aN.jk)]^ 

E[N,{k)] + aE[N,{k)] ' ^'^ 

respectively, where E[Nj{k)] ~ jV^'^Y^jJ^y'Nji^), and f^ji^) is the number of photons 
observed in the detection area si^det,] in the fc-th image. [See Subsection ( 13.41 )1. 

e) Evaluation of the optical losses 

The actual value for the estimated quantum efficiency r\j subsumes also losses due to the 
crystal, the lens and the mirrors. Thus, the value of the quantum efficiency of the CCD is 
obtained as the ratio between the estimated value for the quantum efficiency r\j and the 
transmission on the y-channel, i.e. 

^true.j = —, (8) 

with j — s, i. Tj should be evaluated by means of an independent "classical" transmittance 
measurement. As in this paper we present just a proof of principle of the proposed tech- 
nique, we will not discuss this transmittance measurements anymore in this paper. Thus, 
instead of providing the quantum efficiency of the CCD standalone, the results presented 
in the following can be interpreted as the quantum efficiency of the whole optical system 
before the CCD, including the CCD itself. 

3.1. determination of CS 

The single shot image is stored as a "superpixel" matrix (Fig. |2]i. We select two equal rect- 
angular regions ^/g,,., and sz/^ier.i belonging to the signal and idler portion of the image and 
containing a certain number of "superpixels". In our case, we choose 9 mm^ areas that span a 
wavelength bandwidth of the order of 1 Onm around the degeneracy at 7 1 Qnm. Fixing the region 
■s^det.s, we evaluate the noise reduction factor Ospatiai ^s a spatial average on the pairs of conju- 
gated "superpixel" inside the two regions, in function of the position of the center of the region 
si^det,i of coordinate B, = (^i , ^2)- In this way we obtain a matrix of values of the spatial average 
Osptitiai (^ ) as a function of the center t, of the detection area J^det,i- 
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Fig. 2. Determination of center of symmetry. The right hand side presents a typical image 
obtained by the CCD in the working condition. The images is separated in two portions, 
one collecting the light from the signai beam (H-polarized) and the other the light of the 
idler beam (V-polarized). Two regions Rs and Ri (squares) are selected, and subtracted 
pixel by pixel. The spatial noise of the difference is reported in the inserts (right-bottom) in 
function of the position t, = (^1,1^2) of the idler region. 



The result obtained by the analysis of a typical image are presented in Fig.|2l where a binning 
6x6 of the physical pixel has been applied and the number of photon per superpixel is ~ 1700. 
The presence of correlations around ^ = is represented by a deep in the values of the spatial 
average of Ospatiah whereas far from the minimum the correlation decrease, because conjugated 
pixels no more detect the correlated photons. Thus, this measurement of the spatial correlation 
allows to determine the best position of J^jet.i with respect to s^det.s^ ^nd hence the center of 
symmetry of the correlation. The size of the correlation deep represents the coherence area that 
in our experiment is s^coh ^ 120 x 120(/Xm)^. 

The superpixel size used for the measurement must realize a tradeoff between the visibility 
of the correlation, and the final uncertainty in the CS determination. On one side, as discussed 
in Sec|2] for a good visibility of the quantum correlation (and to increase the signal to electronic 
noise ratio) the "superpixel" area should be much greater than A^„h- On the other side, a small 
pixel size will lead to a small uncertainty. We found that in our setup the best choice is a pixel 
size equal to the coherence area, that can be obtained by a 6x6 binning, as in Fig.|2] 

From the discussion above, it is clear that the measurement described allows a positioning 
of the two correlated (symmetrical) regions within a single "superpixel" uncertainty, while the 
center of symmetry is identify within half a "superpixel". However, it can be demonstrated that 
even a shift of a small fraction of "superpixel" with respect the real CS determine a increasing of 
the noise reduction factor ||2]|9l- In practice, the optimization of the NRF by micro positioning 
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Fig. 3. Noise reduction factor in function of tine detection area, x-axis reports tlie linear 
size of tfie detection areas in ^m. Tine two experimental curves refer to the noise reduction 
factor Oa witfi (red), and without (black) background correction. 



of the CCD allows to determine the physical center of symmetry within a final uncertainty less 
than 1/10 of the "superpixel", hence of the coherence area in our setup |l2]|9l- 

3.2. Determination of the minimum size of s/cif^ j 

In the previous section we showed that CS can be determined with good precision, and the 
regions s^dei.s ^nd s^det.i ^6 consequently fixed to be strictly symmetrical and correlated even 
"locally" (i.e. for each pair of twin spatial modes). 

However, as pointed out in Sec. |2] perfect detection of quantum correlation requires the 
condition ./^spatial = '^det.j/ s^coh ^ 1 . In order to establish when this occurs, we measure the 
noise reduction factor in function of the detection area Oa{Adet), i-e. in function of the number 
of spatial modes detected. The measurement has been performed on a set of ,y^ = 4000 images 
using a binning 12x12 that means a superpixel size of 240 x 240(jLim)^. Here, we define 
A^; [k) the total number of photons detected in the region s^det.j of the ^-th image of the set. 
Oai'S^det.j) has been evaluated according to Eq. (|7]), over the set of 4000 images. Here, the 
single determination of the noise reduction factor is obtained by subtracting the total numbers 
of photons detected in the two large regions s^det,j of the same image, and the statistics is 
obtained over the set of images. The results are reported in Fig. [3] 

As expected, Oa{^detj) is a decreasing function reaching an asymptotic value for si^det.j > 
£/q — 1440 X 1440(/iOT)^ that corresponds roughly to a number of spatial modes larger than 150. 
Therefore we can affirm that, working with detection area larger than ^/q allows to match the 
condition ^spatial = ^det,jl ^coh ^ 1- This prevents possible bias in the estimation of quantum 
efficiency due to the uncertainty in the propagation direction of correlated photons. 

5.5. Experimental contributions to the excess noise 

Other sources of experimental noise lead to systematic issues on the experimental values of 
Ga and consequently on the quantum efficiency obtained by Eq. (|5]l. In the following we will 
address this problems. 



The largest contribution to the excess noise that we have in our setup is related to the in- 
stability of the Q-Switch laser pulse. In particular, we observed a fluctuation of the energy 
pulse-to-pulse of more than 10%. The power P of the pulse is directly related to the mean value 
of photons per mode jx oc sinh'(co«j'f * Vp)- Therefore, the temporal statistics of the PDC emis- 
sion is drastically influenced by the pump power fluctuation. In particular jj. is not a constant. A 
temporal statistics on many pulses (many images) will be characterize by a mean value /T and 
variance V (/i). It can be demonstrated that the contribution of the pump fluctuation modifies 
the expected value of the noise reduction factor with respect to Eq.llUi in the form 






Ar+i + ^(i + .^,„ 



(9) 



From this equation it is clear that the instability of the pump generates a contribution to the 
amount of excess noise that can be relevant, since it include a factor ^tot- U Following the 
same argumentation of Section |2l the effect of this enhanced excess noise can be suppressed 
if the losses on the two beams are a posteriori compensated, by evaluating Ga instead of a. In 
this case da reduces again to Eq. (|5]). 

An other important source of excess noise is the background generated by the electronics 
of the CCD (digitalization, dark counts) and from the straylight, mostly caused by the fluo- 
rescence of the filter and mirrors used for stopping the pump, and residual of the pump itself. 
The first contribution, the electronic one, depend on the level of binning, and can be considered 
independent with respect the thermal noise of the PDC, and straylight noise. In principle, also 
the straylight noise is uncorrected with respect the thermal fluctuation of the PDC light, and 
can be represented by a poissonian-like statistics. However, some correlation between the stray- 
light and the PDC emission is introduced by the fluctuation of the pump, since when the pump 
pulse is more(less) energetic both the PDC and straylight increase (decrease) accordingly. Any- 
way it can be demonstrated that even this correlation cancel out in the difference of the photon 
number when the transmission of signal and idler path are balanced. We define the number of 
counts A'^ ,. registered in the region /?^. ;,■ expressed as the sum of the PDC photons A'^^ /, and the 
background Mj /, . 

The expression linking the quantum efficiency to the expectation value of measurable quan- 
tities in presence of background, the analogous of Eq. (|5]l, is 

_ {5HN',-aBN:))-{5HM.-asM.)) 1 

where Oig = -^ = /jlZiJ) ■ ^^^ background and its statistics can be measured easily and inde- 
pendently, by collecting a set of ^ images when the PDC is turned off, just by a 90" rotation 
of the crystal. Following the same formalism of point d) of the experimental procedure, a^ is 

obtained as; 

^ E[N:{k)]-E[M,,{p)] 

"' E[N',{k)]~Emp)y ^ ^ 

where E[N'j{k)] = ^^^L'ktiN'j{k), and E[Mj{p)] = ^-^Y:f=\Mj{p) represent the experi- 
mental determination of the quantum expectation values (M) and (My) respectively. At the 
same time, aa.B in Eq. dTol ) is obtained by the following experimental estimates: 

{8\N[-aBf^[))^Em{k)~aBf^\{k)f\-E{f^[{k)-aBf^[(k)f (12) 

(52(M, - anM^) ^E[{M,{p) - aBMi{p)f] -E[M,{p) - aMi{p)f 



' The order of magnitude in our experiment is estimable as large as 2r\+^,otV {l^) l~P = V {{N, + Nj)) / {N, +Ni) ~ 
5 ■ 10' (see Tab. [1] first two columns). It is four order of magnitude larger than the excess noise due to the thermal 
fluctuations. 



We also mention that, for each measurement performed in the present work, the images 
affected by cosmic rays have been discarded by a proper algorithm. 

3.4. Efficiency estimation and uncertainty evaluation 

In the previous Subsections we have presented the procedure to obtain an unbiased estimation of 
the detection efficiency Tjj by means of appropriate positioning and sizing the detection regions, 
a posteriori balancing of losses and background contribution analysis. In this subsection we 
present the experimental estimation of 7]^ as well as its uncertainty budget. 

The sizing parameter chosen for this experimental proof of principle are: areas of detection 
■s^detj = 2400 X 3840(^m)^ (5x8 superpixels of size 480 x 480(^m)' obtained by a binning 
24x24 of the physical pixels) corresponding to about 640 x ^0/7 • The estimated value of 77, is 
obtained by the inversion of Eq. ( fTOj i. where ag is calculated according to Eq. (fTTT i and aa,B 
performing the substitutions in Eq.s (fT2] i on the basis of ^ images with PDC light, and ./# of 
background light. 77, can be obviously evaluated as 77, = 0:5775. 

Once we intend to provide an uncertainty associated to the estimated value of 77j (and 77,), we 
should repeat the experiment iF-times, i.e. we should collect jF • ,yK images with PDC light, 
and iF • ^ with background Ught. Thus the estimated value of Ub and Ga.B can be obtained 
from 

a\e.s,i,n = 2f-'Y.a{l) aa\esrin,^2f-'Y.aa{l) (13) 

1=1 1=1 



and the associated uncertainty is obtained following the guidelines of Ref. II32I . Specifically 
the uncertainty propagation is performed on the 2./K + 2^ measured quantities, namely A^^' {k), 
N[{k'), Ms{p), Mi{p') with k,k! = \,...,^ and p,p' = \,...,Ji . We accounted for the correla- 
tions between N[{k) and N[{k') when k~k', due mainly to PDC light, and between Ms{p) and 
Mi{p') when p = /?', due to pulse-to-pulse laser energy instability. Thus, it is assumed that there 
is no correlation between measured quantities in different images. In our experiment jV ~ ^ 
and JT = 500 and jT = 8. 

The estimated values of a^ and Oas together with their uncertainties are presented in Table 1 , 
where, for the sake of comparison, are shown also the corresponding values without background 
subtraction (a and <?« respectively). According to Eq. (fTOl i we obtained 77j = 0.613 ±0.011. 

We underline that the estimated value of 77^ corresponds to the efficiency of the whole quan- 
tum channel including the CCD -not only the efficiency of the CCD itself-, and that the un- 
certainty associated to 77^ accounts only for the statistical contributions due to the fluctuations 
of the A^-s and M-s measured quantities (Type A uncertainty contributions, according to Ref. 

ESI). 

Typically, further non-statistical uncertainties contributions (Type B 1321 ) should be ac- 
counted for in a complete uncertainty budget. For example, in this case a proper evaluation 
of Ub with small uncertainty is mandatory to nullify the effect of the excess noise in Eq. ^, i.e. 
to ensure the validity of Eq. (fTOl i. A Type B uncertainty contribution associated to the nullifica- 
tion of excess noise term should be considered. In Section [373] we showed that the excess noise 
due to the pulse to pulse instability is of the order of 5 • 10^ (see Footnote). From Eq. (|9|l it turns 
out that the condition for neglecting that term is 77^/(477^) • 5 • lO"* <C 1, that in our case means 
77_ ^ 10^^. The uncertainty on a^ presented in table [T] shows that the balancing of the two 
arms can be performed within 10^^, equivalent to the condition 77_ — 10^^. Thus, the possible 
contribution to the uncertainty due to this term is less than 10^^. 

In the determination of Oa.B this uncertainty contribution is several order of magnitude below 
the statistical (Type A) contributions, thus, absolutely negligible in the complete uncertainty 
budget. A larger Type B contribution comes from the bias in the determination of the center of 
symmetry that in our present experiment is 1/10 of the coherence area in both the coordinates 
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Table 1. Experimental estimates (row 2) and their uncertainties (row 3). Column 1-2 show 
the mean counts in the detection area .c/dei,s of the single images and the mean square 
of fluctuation S~N!.{k) = {N!,{k) - E[N'^])-. Column 3-4 report the same values related to 
the background images. Column 5-6 present the experimental values of the compensation 
factor a and the value corrected for background counts. In column 7-9 are reported the 
raw noise reduction factor, the one after compensation of the losses, and finally the one 
after compensation and background correction respectively. 



in the detection plane. Following the argumentation at the end of section |3TT| it generates an un- 
certainty of 1.5%. The actual relative uncertainty in the estimation of Tjj is 2.3%, it is expected 
that it can be easily reduced of more than one order of magnitude increasing the number of 
collected images, as far as the bottleneck of Type B uncertainty contribution is reached. 



4. Discussions and Conclusions 

As a test of consistency for the theoretical model at the basis of the proposed CCD calibration 
technique, and, in particular, for the associated uncertainty model we evaluate the statistical 
uncertainty associated to the mean values of Ub and cJa,B, according to 
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(14) 
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1=1 



obtaining a good agreement with the estimated uncertainty. 

Furthermore, we observe that according to the principle that the accuracy of a measurement 
depends on the measuring time (in our case the number of acquired images), and not depend 
on how the data are arranged, we verified that the final uncertainty on the mean values is not 
influenced by the different possible choices of 3f and ^/K, provided that total number of images 
3f ■ o/K = const. On the contrary, the standard deviation of the populations, is a function of ^, 
as it is shown in Fig. ID 

We note that TablT] shows that (M,) is 5% of the total counts (A^^') although the weigh of 
background correction is the 15% of the estimated noise reduction factor. Nevertheless, the un- 
certainty on the value of Oa.B is just slightly influenced by the background correction. For the 
sake of completeness we observe that the electronic noise contributes with a standard deviation 
A ^ 70 counts. In a forward-looking perspective, for pushing the uncertainty on this measure- 
ments at the level of the best values obtained in the single photon counting regime, it would be 
very important to reduce the straylight. Actually it is possible to design a different experimental 
configuration that limits the fluorescence of the pump. 

In conclusion we have proposed and demonstrated experimentally a new method for the 
absolute calibration of analog detector based on measuring the sub-shot-noise intensity corre- 
lation between two branches of parametric down conversion. The results on the calibration of 
a scientific CCD camera demonstrate the metrological interest of the method, that could find 
various applications, starting from the possibility to give a key element for redefining candela 
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Fig. 4. Standard deviation of Oa.B on a number iF of independent measurements. In each 
measurement a value of OaM is obtained by formula UOI l with the substitutions ( 112b . using 
a set of ^/K images. Although the standard deviation of the population decrease for larger 
^, the uncertainty on the mean value, only depends on the total number of images 2f ■ jV 
[See text for further discussion]. 



unit II33I in terms of photo-counting. 
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